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ABSTRACT 

We present results of long-slit spectroscopy in several slit positions that cover different 
morphological structures of the central parts of three bright Galactic H II regions: M8, 
M17 and NGC 7635. We study the spatial distributions of a large number of nebular 
parameters such as the extinction coefficient, line fluxes, physical conditions and ionic- 
abundances at the maximum spatial resolution attainable with our instrumentation. 
Particularly, our goal is to study the behaviour of the abundance discrepancy factor of 
2+ , ADF(0 2+ ), defined as the logarithmic difference of the 2+ abundances derived 
from collisionally excited and recombination lines. We find that the ADF(0 2+ ) remains 
fairly constant along the slit positions of M8 and M17. In the case of NGC 7635, we 
only detect the O n recombination lines in the integrated spectrum along the whole 
slit, where the ADF(0 2+ ) reaches a remarkably high value of about 0.59 dex. We 
compare our results with previous ones obtained for the Orion Nebula. We find several 
evidences that suggest the presence of a candidate to Herbig-Haro object in M8. 
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1 INTRODUCTION 

The study of the elemental abundances in H n regions is an 
essential tool for our knowledge of the chemical evolution 
of the universe. Traditionally, ionic abundances relative of 
the elements heavier than He have been determined from 
the strong collisionally e xcited lines (CELs). More than 20 
years ago, iFrencbl (| 19831 ) obtained the first determination 
of the C 2+ /H + ratio derived from the faint recombination 
line (RL) C n 4267 A for a planetary nebula (PN), find- 
ing that it was several orders of magnitude larger than the 
abundance obtained from the CELs of th is ion. Later, this 
resul t was confirmed in o ther PNe (e.g. iRola Sz Stasinskal 
1994; Matins & Liu 199 9]). A s imilar qualitative result was 
also found by iPeimbert et ail (jl993al ) for the 2+ /H+ ra- 
tio in the Orion Nebula: the abundances obtained from the 
flux of the faint RLs were higher than those derived using 
the standard method based on CELs. Currently, this obser- 
vational fact is a classical problem in the understanding of 
the physics of photoionized nebulae known as "Abundance 
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Discrepancy" (AD) problem. This disagreement is quanti- 
fied by means of the Abundance Discrepancy Factor (ADF) 
that can be defined as the ratio, or the logarithmic differ- 
ence, between abundances of a same ion derived from RLs 
and CELs. In the case of the 2+ /H + ratio, the ADF has 
rather similar values between 0.1 and 0.3 dex for extragalac- 
tic and Galactic H n reg i ons ( e .a. iGarcfa-Roias et al.l 120051 : 
iGarcfa-Roias fc Estebanl 120071 : lEsteban et alj|2009h . while 



for PNe the ADF shows a much wider range of values, be- 
coming substantially larger in some object (e.g. iLiu et all 
1200(1 120061 : iTsamis et al.ll2004 120081 ). 

What causes the AD problem is nowadays debated. On 
the one hand, the predictio ns of the temper ature fluctua- 
tions paradigm proposed by iPeimbertl l| 19671 ). and charac- 
terized by the mean square of the spatial distribution of 
temperature -the so-called temperature fluctuations param- 
eter, t 2 - seems to e xplain the ADF observed i n H 11 regions, 
as it is argued by iGarci'a- Roias fc Estebanl (|2007D . Under 
this scheme, the AD problem is a direct consequence of the 
different temperature dependence of the emissivities of the 
lines used: in the case of CELs, it depends exponentially 
on the electron temperature, T e , of the ionized gas, while 
the emissivity of RLs has a power law temperature depen- 
dence, similar to those of the Balmer lines used as refer- 
ence to determine the ionic abundance ratio relative to H + . 
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On th e other hand, the hypothesis suggested by iLiu et al.l 
(2000), where most of the emission of RLs come from a 
cold hydrogen-poor component immersed in the ambient gas 
which emits the bulk of CELs, seems to solve the AD prob- 
lem in PNe with high ADF values. This hypothesis is based 
on the observed fact that certain PNe contain well resolved 
H-deficie nt knots which are strong metall ic RL emitters {e.g. 
Abell 30. lHarrington fe Feibelmanlll984l ). 

The existence and origin of the temperature fluctua- 
tions is controversial because high values of the t 2 parame- 
ter are not reproduced by standard photoionization models 
(|Kingdon fc Ferland|[l995l ; iRodrfguez fc Garcfa-Roiasll201(]| ) 
and additional mechanisms are proposed in order to explain 
the presence of temperature fluctuations (see revisions of 
lEstebanll2002l ; iPeimbert fc Peimbertl2006l ) . In the same way, 
new scenarios try to find a solution to the AD problem in H II 
regions putting forward new physical mechanisms. This is 
the ca se of t he hypothesis pr e sented bv lTsamis fc Pequignod 
l|2005l) and IStasiriska et~afl (|2007t ). Based on the chemi- 
cal m odel for the heavy-elements mixing of iTenorio-Taglei 
(1996), these authors proposed the presence of two compo- 
nents of different chemical composition and physical condi- 
tions in H n regions. The component responsible of most 
of the emission of RLs consists of cold metal-rich droplets 
from supernova ejecta still not mixed with the ambient gas 
of the H II region where most of the CEL emission would be 
produced. Then, note that assuming a chemically inhomoge- 
neous model in H II regions in order to explain the AD prob- 
lem, the abundance derived from RLs and CELs would be 
upper and lower limits, respectively, o f the real abundance of 
the ionized gas (Stasiriska ct al. 2007). Recen tly, a new sce- 
nario have been proposed by lErcoland (|2009h based on the 
existence of high-density quasi-neutral clumps -embedded 
in the nebular gas ionized by the extreme-ultraviolet (EUV) 
radiation- which are ionized mainly by the X-ray emission 
from the central star. Under this scheme, the CEL emission 
mainly comes from the region ionized by the EUV radiation 
(E region), while the RLs are emitted in different propor- 
tions from the clumps (X region) and the E region. In this 
sense, the abundances of the E region would be representa- 
tive of the nebula and those from CELs would be easier to 
correct than RL one s . Con trary to the model propose d by 
iTsamis fc Pequignot] (12005 ) andlStasinska et all (|2007h . the 
nebular model of lErcolano (|2009l ) has homogeneous abun- 
dances. 



In a previous paper |Mesa-Delgado et al.l 1200*81 ) we ex- 
plored the behaviour of the AD at small spatial scales and 
its dependence on different nebular parameters and physical 
conditions in the Orion Nebula in order to shed light on the 
origin of the AD problem. In that study, we used long-slit 
spectroscopy at spatial scales of l'.'2, finding high ADF(0 2+ ) 
values related to the presence of Herbig-Haro (HH) objects 
and temperature spikes at the position of the protoplanetary 
discs (proplyds). A subsequent detailed analysi s of HH 202 
was carried out by iMesa-Delgado et al. (2009a) using inte- 
gral field spectroscopy confirming a high ADF(0 2+ ) value 
at the main knot of HH 202 in agreement with the results 
of Mesa-Delgado et alj ( 2008). Another important result of 
IMesa-Delgado et al l l|2009al ) was the obtention -for the first 
time in an H II region- of a map of the Balmer temperature 
and of the temperature fluctuations in the observed field, 



finding no correlation between the ADF(0 2+ ) and the t 2 
parameter. 

Following the same goals and methodology of 
M esa-Delgado et all (|2008l ). in this paper we have used long- 
slit spectroscopy at intermediate spectral resolution in order 
to study the spatial distribution of the ADF(0 2+ ) and other 
main nebular properties as well as their relation with the 
local morphological structures (e.g. density condensations, 
ionization fronts or HH objects) in other bright Galactic H II 
regions, namely, M8, M17 and NGC 7635. 

After the Orion Nebula, M8 and M17 are probably the 
most studied Galactic H II regions. M8 forms a blister of 
photoionized material on the surface of a giant molecular 
cloud. Near the optical center the region with the high- 
est surface brightness of the nebula is found, the Hourglass 
(HG) region. This region is mainly ionized by the O star 
Herschel 36 (Her 36), while t he stars HP 165052 and 9 Sgr 
ionize the rest of the nebula (jWoodward et al.| [T986). M17 
is a cavity with a V shape and the V opening in the line of 
sight. The main ionization source of M17 is a group of 03-04 
stars, which belong to the open cluster NGC 6618, loca ted 
in the dark bay of the nebula l|Hanson fc Contil 1 19951 ). A 
singular characteristic of M17 is its rather high ionization 
degree -0 2+ /0+ ratio- in comparasion with other Galac- 
tic H II regions. The chemical composition of M8 and M17 
have been widely studied by several authors in al l spec- 
tral ranges from low to high spectral resolution (e.g. iRubinl 
" Peimbert fc Coster o 1969; Sanchez fc Peimbertlll99'll; 



1969; 



Rodri guez 1999; Esteba n et all Il998l ; ITsamis et all |2003) 



Base d on high resolution and d eep echelle spectrophotom- 
etry, iGarcfa-Roias et al.l (|2007l ) have provided a complete 
revision of the chemical abundances of these H II regions 
using CELs and RLs. Nevertheless, we have not found in 
the literature detailed studies about the spatial behaviour 
of the nebular pr o pertie s of these regions, excluding that 
of IPeimbert et all \ 19921 ) along of 17 areas of Ml 7 at low 
spectral resolution. Our third region, NGC 7635, is not a 
classical H II region. This nebula is an interstellar bub- 
ble formed by the interaction of the stellar wind of the 
06.5 Illf star BD+60 2522 with the surrounding interstellar 
medium. Recently works assume that the ram pressure of 
the stellar wind is balanced by the surrounding gas pres- 
sure due to the similarity between the velocities of the 
the molecul ar cloud and the bright nebulosity places o f 
the nebula (|Christopoulou et al.l [l995l ; iMoore et al.ll2002l ). 
The study of physical conditions and chemical abundances 
in NGC 7635 have been restricted to s ome selected zone s 
dTalent fc Dufourlll979l ; lRodriguedfl999l ; IMoore et al.ll2002l ) . 
IMoore et al. (2002) have been the first in exploring the spa- 
tial distributions of several bright emission lines -[O ill] 5007 
A, Ha and [N n] 6584 A- along the set of knots located at 
northwest of the central star and the rim of the bubble. 
These authors also obtained the first density spatial profile 
along the slit position that covered the knots. 

In Section [2] we describe the observations of the Galac- 
tic H II regions, the reduction procedure and the extraction 
of the one-dimensional spectra. In Section [3] we enumerate 
the selected emission lines and describe the procedure used 
to measure the fluxes and the extinction correction applied 
for each nebula. In Section [3] we describe the method used 
to determine the physical conditions and the ionic abun- 
dances from both kinds of lines, CELs and RLs. In Section[5] 
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Figure 1. Observed slit positions over the central parts of the three Galact ic H II regions. The s lits only show the total extraction area, 
while the total slit length is 3.'8. The positions of the candidates to proplyd llStecklum et al.lll998l) and HH object (see i|7.3l l are indicated 
as UC and HH?, respectively. The images are combinations of exposures taken with different narrow-band filters. For all regions, emission 
from [O III] is shown in blue, emission from Ho is shown in g reen and emission from [S II] is shown in red. M8 and NGC 7635 images 
are sections of combinations of WFPC2 images - tCauld jl997T) and lMoore et ah | J2002I) , respectively. The M17 image is a section of the 
original mosaic obtained by the amateur astronomer I. de la Cueva Torregrosa. 



we present and discuss the spatial distributions along the 
slit positions of several nebular quantities for each nebula. 
In Section [6] we show the physical conditions and the ionic 
abundances for the individual extractions of NGC 7635 as 
well as discuss the puzzling abundance pattern found in this 
object. In Section [7] we compare the results on the light of 
the different lineal spatial resol ution attained in this paper 
with that of the observations of Mesa-Delgado ct al . (2008) 
in the case of the Orion Nebula. We also present several ar- 
guments for the presence of a new HH object in M8 and 
discuss the possible causes of the high ADF(0 2+ ) found in 
NGC 7635. Finally, in Section [8] we summarize the main 
conclusions of the paper. 



2 OBSERVATIONS, DATA REDUCTION AND 
EXTRACTION OF THE 
ONE-DIMENSIONAL SPECTRA 

Long-slit spectra at intermediate spectral resolution were 
obtained on 18 July 2007 and 11 May 2008 using the Inter- 
mediate dispersion Spectrograph and Imaging System (ISIS) 
at the 4.2m William Herschel Telescope (WHT) in the Ob- 
servatorio del Roque de los Muchachos (La Palma, Spain). 
Two different CCDs were used at the blue and red arms 
of the spectrograph: an EEV12 CCD with a configuration 
4096x2048 pixels with a size of 13.5/xm per pixel in the blue 
arm and a REDPLUS CCD with 4096x2048 pixels with a 
pixel size of 15/xm in the red arm. The spatial scale was fy.'20 
pixel -1 and 0'.'22 pixel -1 in the blue and red arm, respec- 
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tively. The slit length was 3'8 and the slit width was fixed 
to 0'.'98. The R1200B grating was used in the blue arm and 
the R316R one in the red arm. These gratings gave an effec- 
tive spectral resolution of 0.86 and 3.81 A for the blue and 
red arms, respectively. The blue spectra covered the spec- 
tral range from 4220 to 5080 A and the red one from 5320 
to 8100 A. The observation nights were photometric and the 
seeing during both observations was between of 0'.'5 and 0'.'8. 

The Galactic H II region M8 and the wind-blown nebula 
NGC 7635 were observed on 18 July 2007. In the case of M8, 
we observed two slit positions (POS1 and POS2) centered 
at the HG region with different position angles (see Fig. (TJ 
the HG region is indicated by a grey section over the slit 
of POS1). These positions were chosen in order to cover 
different morphological struct ures such as a candidate to 
proplyd (|Stecklum et al.lll998l ) marked as UC in Fig. Q]and 
the prominent HH 870. Unfortunately, during the analysis 
of these data, we noticed that probably the proplyd was not 
covered by the slit as we expected due to the difficulty of 
possitioning the slit so close to the bright star Her 36. In 
the case of NGC 7635, a single slit position was used to 
cover the bright knots Kl and K2 as well as the rim of the 
bubble (see Fig. [T]). On 11 May 2008 a single slit position 
was observed over the Galactic H n region M17 covering a 
high surface brightness zone. For all objects, large and short 
exposures were taken at each slit position and spectral range 
in order to achieve a good signal-to-noise ratio in the faint 
C n and O II RLs and to avoid saturation of the brightest 
emission lines. The journal of observations can be found in 
Table [1] where we present the coordinates (RA, DEC) and 
the position angle (PA) of each slit position observed as well 
as the total exposure times. 

All CCD frames were reduced using the standard irai{]] 
twodspec reduction package to perform bias correction, 
flat-fielding, cosmic-ray rejection, wavelength and flux cal- 
ibration. The wavelength calibration was carried out with 
a CuNe+CuAr lamp. The absolute flux calibration was 
achieved by observations of the standard stars BD+33 2642, 
BD+28 4211 and BD+25 4655 for the first night and 
BD+25 3941, Feige 34 and BD+33 2642 for the second one. 
The error of the flux calibration is of the order of 5 per cent. 

The extraction of the one-dimensional spectra was done 
using an iraf script based on apall task , follo wing the 
procedure explained in iMesa-Delgado et al.l l|200ct ). We ex- 
tracted the apertures for each region applying the bidimen- 
sional fit to the spectra of a standard star (with positional 
coordinates more similars to the object) used for flux cali- 
bration on the bidimensional image. In all cases, we adjusted 
a third-order Chebyshev polynomial obtaining a typical rms 
between 0.05 and 0.08 pixels. The slit center in the red arm 
was some pixels displaced with respect to the slit center in 
the blue arm; this effect was also corrected in the extrac- 
tion procedure, and later verified from the alignment of Ha 
and H/3 spatial profiles, ensuring the same spatial coverage 
in both ranges. We also discarded apertures located at the 
edges of the CCD. 



1 IRAF is distributed by National Optical Astronomical Obser- 
vatories, operated by the Associated Universities for Research in 
Astronomy, under cooperative agreement with the National Sci- 
ence Foundation 



For M8 we extracted apertures of 6 pixels in the blue 
arm, and 5.45 pixels in the red one, in the spatial direction, 
which corresponds to an effective spatial resolution of l'/2. 
This size was chosen as a compromise to have the maximum 
attainable spatial resolution (only slightly larger than the 
mean seeing of the night) and a good signal-to-noise ratio 
in the auroral lines of [O ill] and [N n]. However, larger 
extractions of 4'.' 8 angular size were necessary to achieve a 
good flux measurement in the fainter O II RLs. Then, for the 
slit position 1 and considering the discarded apertures, we 
have obtained a final number of 150 apertures with a spatial 
resolution of 1'.' 2x0'.' 98 and 38 apertures with a resolution 
of 4'.'8x0'.'98 for the measurements of O II RLs. Similarly, a 
total number of 146 and 37 apertures were extracted, respec- 
tively, for the slit position 2. The smaller number of aper- 
tures obtained in this slit position is due to the fact that we 
discarded 4 apertures which were severely contaminated by 
stellar emission from Her 36. 

In the case of NGC 7635 slit position, we extracted aper- 
tures with an angular size of 3". A total of 30 apertures with 
an area of 3"x0'.'98 were used. We also extracted individual 
integrated spectra of knots Kl and K2, and the rim of the 
bubble with areas of 12'.'4x0'.'98, 8'.'5x0'.'98 and 10'.'4x0'.'98, 
respectively. We notice the detection of the C II 4267 A RL 
in the three extractions. 

For M17, apertures extracted with an angular size of 
1'.'2 was sufficient to obtain a high signal-to-noise ratio in 
the auroral lines and RLs. We have obtained a total num- 
ber of 185 apertures extracted with an individual area of 
l'.'2x0'.'98. 

Additionally, for each slit position, we extracted a one- 
dimensional spectra collapsing the whole slit -the sum of all 
the individual apertures. These are designated as "whole 
slit" spectra. The area covered by the "whole slit", and 
therefore the total extraction area, is shown in Fig. [T] 



3 EMISSION LINE MEASUREMENTS AND 
REDDENING CORRECTION 

The emission lines considered in our analysis were selected 
according to the following criteria: 

• Hi lines - Ha, H/? and H7-, which are used to compute 
the reddening correction and to re-scale the line flux ratios 
of the red spectral range with respect to the blue one. 

• CELs of several species in order to compute the phys- 
ical conditions such as the auroral lines [O 111] 4363 A and 
[N n] 5755 A used to derive the electron temperatures, [S 11] 
6717, 6731 A and [CI III] 5718, 5738 A used to calculate the 
electron density. 

• Other CELs needed to derive different ionic abundances 
(N+, 0+, 2+ , S+ S 2+ , Cl 2+ and Ar 2 + ) 

• Faint RLs of C 11 and O II, which are used to derive the 
C 2+ and 2+ abundances and to compute the abundance 
discrepancy factor, ADF, for 2+ (via a comparison with 
the 2+ abundances derived from CELs). 

Line fluxes were measured applying a single or a multi- 
ple Gaussian profile fit procedure over a local continuum. All 
these measurements were made with the SPLOT routine of 
the iraf package and using our own scripts to automatize 
the process. Due to the local variations of the continuum 
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Table 1. Journal of observations. 











Exposure time (s) 


Target 


RA a 


DEC 


PA (°) 


Blue arm 


Red arm 


M8 POS1 


18 h 03 m 41?40 


-24°22'42'.'7 


3 


7x1200,60 


36x180,60 


M8 POS2 


18 h 03 m 40 a 42 


-24°22'42'.'7 





7x1200,60 


36x180,60 


M17 


18 ft 20 m 42?98 


-16°10'02'.'4 


120 


5x1800,60 


15x500,60 


NGC 7635 


23' l 20 m 44?36 


+61°11'56'.'0 


67.5 


5x1200,60 


8x600,60 



a Coordinates of the slit center (J2000.0). 



around O n RLs and their faintness, these lines were mea- 
sured manually; 

Following Mosa-Dclgado et al. (2008), to accurately 
compute the line fluxes we need to define the adjacent con- 
tinuum of each line using the SPLOT routine. For each se- 
lected line, we define two small spectral zones at each side of 
the line as close as possible and free of any spectral feature. 
Then, the routine fits the continuum between both zones 
and the line profile to obtain the flux. The observational er- 
rors associated with th e flux measurements were e stimated 
following the criteria of iMesa-Delgado et al] (120081 ). The fi- 
nal error of a line flux was computed as the quadratic sum 
of the error in its flux measurement and the error in the flux 
calibration. All line fluxes for a given aperture were normal- 
ized to an H I line of reference, H/3 and Ha for the blue and 
red range, respectively. 

The observed fluxes with respect to their H I line of 
reference were dereddened using the usual relation, 

7(A) _ -F(A) 1n c(Hg)[j(A)-/(y ef )] m 
I(Xr Bf ) F(\ ref ) ' 1 } 

where the reddening coefficient, c(H/3), represents the 
amount of interstellar extinction, /(A) the adopted extinc- 
tion curve normalized to /(H/3) = 0, and A re / the H I line 
of reference. The reddening coefficient was determined from 
the comparison of the observed flux ratio of H7 and Ha with 
respect to H/3 and t he ca se B theoretical ones computed by 
IStorev fc Hummerl (1995) for the physical conditions of T e 
= 10000 K and n e = 1000 cm" 3 . The final c(H/3) was the 
weighted average of the values obtained from each line. In 
appendix [X] we present the dereddened emission line ratios 
and their associatted errors of the main emission lines per 
slit position as well as the observed H/3 flux and the mean 
c(H/3) coefficient (just in the online version). 

In all regions we have a ssumed the extinction law de- 
rived bv lCardelli eTaT] (| 19891 ), which is parametrized by the 
ratio of total to selective extinction, Rv = Av / E(B — V). In 
the case of M17, NGC 7635 and the position 2 of M8, we have 
used the typical value in the diffuse interstellar medium, 
Rv = 3.1. For the position 1 of M8 we have adopted differ- 
ent values of Rv depending on the zones covered by the slit. 
It is we ll known that the mai n ionization source of HG is 
Her 36 l|Woodward et al.lll986l ) and that it shows a consid- 
erably higher extinction than the other zones of M8, which 
are mainly ionized by 9 Sgr and most of their reddening is 
due to foreground interstel lar dust. Consequently, following 
ISanchez fc PeimbertJ l|l99ll ) we have adopted a Rv = 5.0 to 
correct the apertures of position 1 that cover the HG region 
(the area indicated by a darker grey band in the slit position 
1 of M8 shown in Fig. QJ and the typical value Rv — 3.1 in 
the other zones. The use of a higher total to selective extinc- 



tion produces significant changes in the c(H/3) values, while 
the line flux ratios remain almost unaffected by deviations 
from the classical extinction law. In the case of the "whole 
slit" spectra, we have also assumed the typical Rv value 
for the interstellar medium. Finally, in order to produce a 
final homogeneous set of dereddened flux ratios, all of them 
were re-scaled to H/3. The re-scaling factor used in the red 
spectra was the theoretical Ha/H/3 ratio for the physical 
conditions of T e = 10000 K and n e = 1000 cm" 3 . The final 
error associated with the line dereddened fluxes include the 
uncertainties in the flux measurements, flux calibration and 
the error propagation in the reddening coefficient. 

In the different panels of Fig. [5] we present the spatial 
profiles of the c(H/3) determinations as a function of the 
positional measurement along the slit in arcseconds. Here- 
inafter, each positional measurement represents the position 
on the slit of a given quantity obtained for an individual ex- 
traction with origin in the south edge of the original slit in 
the case of M8, in the northwest edge for M17 and in the 
southwest edge for NGC 7635. In all cases, we have found a 
good agreement between our determinations and those avail- 
able in the literature, considering that the measurements do 
not correspond exactly to the same spatial zone and the ar- 
eas covered by the slits are also different. On the one hand, 
Fig s I2I a ) and[^b) show the spatial distribution of c(H/3) for 
slit positions 1 and 2 of M8, respectively, where we can no- 
tice higher c(H/3) values associated with the HG region and 
the zones near Her 36. These values are in agreement with 
pr evious c(H/3) determinations carried out in the HG r egion 
bylSanchez fc Peimbertl (Il99ll ). lEsteban et all (|l999bl) and 
iGarcfa-Roias et al.l (2007), who found values between 0.85 
and 1.0 dex. On the other hand, in Fig. [2jc) we present the 
spatial distribution of c(H/3) along the slit position of M17. 
W e have compared our de terminations with those performed 
bv lPeimbert et alj (l992), whose slit positions 3 and 13 ap- 
proximately coincide with the positional me asurements 60" 
and 1 25", respectively. For those positions. iPeimbert et all 
1 19921 ) obtained values of about 1.72 and 1.45 dex, which 
are in agreement with our determinations considering the 
uncertainties in the c(H/3) determination. Fig.[2jd) presents 
the c(H/3) spatial profile of NGC 7635. In this panel, we 
can see higher c(H/3) values related to the knots and the 
rim. We have com pared our values with pr evio us determina- 
tions obtained bv lTalent fc Dufourl (|l979l ) and I Moore et al] 
(|2002|) at differen t part s of the nebula. Positions 3 and 5 of 
iTalent fc Dufourl (| 19791 ) coincide with th e positions of the 
knot Kl and the rim indicated in Fig.l2Td). lTalent fc Dufourl 
l|l979l ) measured c(H/3) values at those positions of about 
1.21 and 1.07 dex, respectively, which ar e in agreement with 
our determinations. On the other hand. [Moore et all (2002) 
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Figure 2. Spatial profiles of the reddening coefficient, c(H/3), along the slit positions: a) M8 POS1, b) M8 POS2, c) M17 and d) 
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line represents the c(H/3) value obtained for the "whole slit" spectra. The typical error bar is also included. 



obtained a similar value, 1.25 dex, from their slit position 
that covers the knots Kl and those located to the south of 
Kl, but measured a value of about 1.46 dex for the rim of 
th e bubble, which is higher than our determination and that 
of lTalent fc Dufouij <|l979f ) at the same position. The typical 
error in the c(H/3) coefficient presented in Fig. [2] is of about 
0.1 dex for all regions. 

The results shown in Fig. [2] indicate that an important 
fraction of the extinction should come from dust located in- 
side the objects, and that the distribution of such absorbing 
material is not homogeneous. 



4 PHYSICAL CONDITIONS AND CHEMICAL 
ABUNDANCES 

4.1 Physical conditions 

We have determined the physical conditions -electron den- 
sities and temperatures- from the usual CEL ratios and 
using the iraf task temden of the nebular package 
(IShaw fc Dufouil 1 19951) with updated a tomic data (see 
et al.l 120001 ; iGarcfa- Roias et all 120051 ). We have com- 
puted the electron density, n e , from the [S n] 6717/6731 
line ratio and the electron temperatures, T e , from the 
nebular to auroral [O ill] (4959+5007) /4363 and [N n] 
(6548+6584) /5755 line ratios. Although we detect the [CI ill] 



doublet in several apertures, we do not use these lines to ob- 
tain density values due to their large associated errors. The 
spatial distributions of the physical conditions are presented 
and discussed for each region in Sj5] 

Following the same methodology as 

M esa-Delgado et all (|2008l ) for the determination of 
the physical conditions, a representative initial T e — 10000 
K is assumed in order to derive a first approximation of 
n e ([S 11]) -hereinafter n e . Then, we calculate T e ([0 in]) 
and T e ([N 11]), and iterate until convergence to compute 
the finally adopted values using T e ([N 11]) in the density 
calculations. The errors in the physical conditions were 
computed by error propagation on the analytical expres- 
sio ns of n e by Castaficda et al.| 1 19921 ) and those of T e given 
by lOsterbrock fc Ferlandl ( 2006 ) (their e quations 5.4 and 
5.5). A lthough the expression derived by ICastaneda et al.l 
is only valid to a limited range of densities lower 
than 10 4 cm -3 , an d uses the old a tomic data for S + from 
the compilation by iMendozal (Il983l ). it seems adequate for 
an estimation of the errors in the physical conditions. 



4.2 Ionic abundances from CELs and RLs 

We used the iraf package nebular in order to derive ionic 
abundances of N+, + , 2+ , S+, S 2+ , Cl 2+ and Ar 2+ from 
CELs. We have assumed no temperature fluctuations in the 
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Figure 3. Sections of the spectrum around the emission lines of multiplct 1 of O II. Panels a) and b) correspond to apertures of 4" 8 wide 
extracted from slit positions 1 and 2 of M8, respectively, extending from positional measurements 137'.'2 to 142". Panel c) corresponds to 
an 1"2 wide aperture extracted from the slit position of M17, which covers from 153'.'4 to 154'.'6. Panel d) shows a section of the "whole 
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ionized gas (t 2 — 0) and a two-zone scheme, adopting the 
T e ([N n]) to derive the abundances of singly ionized species 
and T e ([0 ill]) in the case of doubly ionized species. The elec- 
tron density obtained from the[S n] line ratio was adopted 
for all ionic species. The errors in the ionic abundance deter- 
minations were calculated as the quadratic sum of the inde- 
pendent contributions of temperature, density, and line flux 
uncertainties. The spatial distributions of the ionic abun- 
dances of some species are presented and discussed for each 
region in ij5] 

The wavelength range covered with the blue arm did 
not include the bright [O n] 3726, 3729 A lines. In the case 
of NGC 7635, it was necessary to determine the + abun- 
dance (see [HU; to do that, we performed a proper subtrac- 
tion of the telluric emission in the [O n] 7320, 7330 A lines. 
The subtraction was achieved using re-scaled pure telluric 
measurements from sky extractions of different sizes of long 
exposure spectra of a compact H II region observed dur- 
in g the same night (the data of this object were published 
in iMartm-Hernandez et all 120081 ). This procedure adds an 
uncertainty between 6% and 8% higher in the flux measure- 
ments of [O n] lines. On the other hand, in the cases of M8 
and M17 we consider not necessary to determine the + 
abundance for our aims. 

We have detected and measured pure RLs of the mul- 
tiplet 1 of O II (see Fig.© and C II 4267 A. The abundance 
of a heavy element X in the ionization state i + 1, that emits 



a RL at wavelength A is given by 

N(K +1 ) _ A (A) a eff (}ip) /(A) 



AT(H+) 4861 a e //(A) I(H/3) ' ^ 

where I(A)/7(H/3) is the dereddened flux ratio, a e //(A) and 
a e //(H/3) are the effective recombination coefficients for the 
RL and H/3, respectively. Due to the similar temperature de- 
pendence of the emissivities of RLs, the a e ff(Hl3)/a e ff(\) 
ratio is almost independent of the adopted temperatures. 
In our case, under the two-zone scheme, we have assumed 
T e ([0 ill]) to calculate the C 2+ and O 24 " abundances. 

We have detected the C II line at 4267 A in most of 
the apertures extracted at l'.'2 in M17 and M8, and in the 
integrated spectra for the "whole slit", the knots and the 
rim of the bubble of NGC 7635 (see C 2+ abundances 
from RLs h a ve be en calculated using the atomic data of 
iDavev et al] (|2000T ). 

The 2+ abundances were derived when at least four 
lines of multiplet 1 were measured in a given one-dimensional 
spectrum usin g the effective recombination coefficients from 
Storevl (Il994h and the same method that is detailed in 



Esteban et al.l (|l998l ). We have determined the 2+ abun- 



dances correcting for the departure from local thermodi- 
namyc equilibrium (LTE) of the upper levels of the tran- 
sitions of multiplet 1 of O II for densities lower than 
10000 cm" 3 using th e empi rical formulation proposed by 
iPeimbert fc Peimbertl {2005). Abundances determined as- 
suming LTE for the population of the levels do not differ 
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more than a 2% from the NLTE ones because we use several 
of the brightest lines of the multiplet. We have computed 
the 2+ /H + ratio for 89% of the apertures extracted at l'.'2 
in M17, for all the apertures extracted at 4'.'8 in both slit po- 
sitions of M8, and in the "whole slit" spectrum of NGC 7635 
(see p. 



5 SPATIAL PROFILES ALONG THE SLIT 
POSITIONS 

In the following sections we present the spatial profiles of 
several nebular properties along the slit positions of each 
H II region. The selected parameters were: n e , T e ([N n]), 
T e ([0 in]), the flux of several selected lines ([Fe m] 4658 A, 
C II 4267 A, O II 4649 A and [O ill] 4959 A), and the 2+ 
abundances obtained from CELs and RLs. 

5.1 M8 

In Figs U and [S]we show the spatial variations of the nebular 
properties for the slit positions 1 and 2 of M8, respectively. 

The spatial profiles of n e along the slits show a wide 
range of variation with a maximum associated with the HG 
region in position 1 (see Fig. g^i) and the apertures near 
Her 36 in the case of position 2 (see Fig. [5^). The den- 
sity peak in Fig. HJa) associated with HH 870 reaches a 
value higher than 800 cm" 3 . Following the nomenclature 
proposed by lArias et al.l l|2006h . that peak would be more 
specifically related to the knot B of this object. At the HG 
region we find densities between 2600 and 3500 cm -3 , val- 
ues which are in agreemen t wit h previous determinatio ns by 
iGarcfa-Roias et all < |2007h and lPeimbert et al.l (1199301 ). who 
derive densities between 1800 and 3700 cm" 3 . We have no- 
ticed a slight localized enhancement of about 300 cm" 3 in 
the density profile of position 1 around the positional mea- 
surement 155'.'2 (indicated as HH? in Figs [T] and [4]). This fea- 
ture coincides in Fig. [1] with a relatively bright rim crossed 
by the slit and located about 45" to the north of the centre 
of the HG region. In i]7.3l we argue that this feature may be 
a candidate to HH object. In Fig. Oa), the highest density 
value is related to an ionization front located at 15'.' 3 S and 
l'.'5 E from Her 36. In this slit position, we have also de- 
tected local maxima at the positional measurements 65'.'2, 
79'.'6, 151'.'6 and 192? 4. On the one hand, the first two max- 
ima coincide in Fig. [T] with gaseous arcs that surround faint 
stars located approximately at 34" and 49" to the south of 
Her 36. On the other hand, the other maxima are related to 
high density filaments located approximately at 40" and 80" 
to the north of Her 36. The average value of both slit posi- 
tion 1 and 2 obtained from the "whole slit" spectrum, and 
included in panel (a) of Figs |4][5] with dashed line, amounts 
to 1140±220 and 750±150 cm" 3 , respectively. 

The spatial profiles of T e ([0 ill]) and T e ([N n]) are pre- 
sented in the panels (b) of Figs|4]and[5]for the slit positions 1 
and 2, respectively. We do not find relevant features in these 
profiles, excluding perhaps very slight enhancement of both 
temperature indicators at the HG region. However, the tem- 
perature rise at these positional measurements are of the or- 
der of the error bar. We also present the T e ([N Il])/T e ([0 ill]) 
ratio in the bottom panels of Figs[4jb) and[5jb), which is 
basically constant along the slits showing spatial variations 



of the order of the observational errors. In general, T e ([N n]) 
is slightly higher than T e ([0 ill]), a typical result obtained 
for H II regions as a consequence of the harden ing of the 
radia tion field in the low ionization zones (e.g. IStasiriskal 
ll98d V 

One of the main spectral properties of HH objects is 
their strong e mission in [Fe hi] lines due to the destruction 
of dust grains (Mesa-Delgado et al. 2009b) or changes in the 
ionization conditions (jBlagrave et alj|2006l ). In Figs|4jc) and 
[5jc) we have plotted the spatial profile of the dereddened 
flux of [Fe ill] 4658 A with respect to H/3 for the slit po- 
sitions 1 and 2, respectively. In the slit position 1 we can 
see a clear rise of the flux at the position of HH 870, which 
increases by a factor of ~2 with respect to the adjacent 
background gas. This factor seems to extend along the HG 
region where the ionization degree decreases and the elec- 
tron density increases. As in the density spatial profile of 
slit position 1, we have found in Fig. [4fc) an enhancement 
of the [Fe ill] emission line at the positional measurement 
155" 2 where a possible candidate to HH object is located 
(see i]7.3[) . In the case of the slit position 2, we have found 
no evidence of important localized increase in the [Fe ill] 
emission, only a decrease of a factor of 2 around Her 36. 

In Figs QId)-[5f d) an d I2e)-[5le) the spatial profiles of 
the observed flux of C II 4267 A, O II 4649 A, and [O ill] 4959 
A lines along the slit positions are shown, as well as their 
ratio in the bottom panels of the same figures. It should be 
reminded that the extractions of the O II line are 4'.'8 wide. 
Then, each extracted aperture represents the average value 
of four apertures extracted with a spatial size of l'.'2. In or- 
der to perform a proper comparison, we have also plotted 
the observed fluxes of C II and [O ill] lines with extrac- 
tions 4'.'8 wide. The spatial distributions of the pure RLs 
are quite similar in each slit position. However, we can note 
a slight rise of the F(C II 4267 A)/F(0 II 4649 A) ratio at the 
HG and HH 870 regions (see Fig. 2J1). We have observed a 
similar enhancement of the ratio of emission lines of species 
with similar ionization potential to C 2+ (24.4 eV), like Ar 2+ 
(27.6 eV) or Cl 2+ (23.8 eV), with respect to [O ill] 4959 A 
and O II 4649 A. In order to explore the reason of this en- 
hancement, we have determined the spatial distribution of 
the ionization degree which presents an inverse behaviour 
to the C n/O II ratio. Though we have not subtracted the 
telluric emissions from [O n] 7320, 7330 A lines in M8, we 
have observed that the 2+ /0 + decreases between the posi- 
tional meassurements at 70" and 130", where the HG region 
is located, reaching values of about —0.4 dex lower than in 
the rest of the slit position. Therefore, this indicates that 
the variation in the C n/O II ratio can be simply due to the 
decrease of the ionization degree in that zone. On the other 
hand, the point-to-point comparison of the spatial profiles 
of [O ill] 4959 A and O II 4649 A (see Figs gH^) does not 
show clear tendencies. 

The [O ill] spatial distribution of the slit position 1 
and 2 of M8 suggests that the gas covered by both slits 
might not be ionized by the same source. Slit position 1 
shows a [O ill] emission distribution rather similar to the 
density distribution, but position 2 does not. In fact, the 
[O ill] emission shows a peak at the ionization front at the 
south of Her 36 (positional measurement ~95") and remains 
almost constant toward the north of the slit. This area co- 
incides with the dark zone around Her 36. It seems that 
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Figure 4. Spatial profiles of several nebular parameters along the slit position 1 of M8. Positional measurements along the slit go from 
south to north (see Fig. [TJ- The position of the Hourglass (HG) region and the HH object 870 are indicated as well as the typical error 
of some variables. The possible candidate to HH object (see SI7.3H is also marked as HH? in panels (a) and (c). The horizontal dashed 
line in some panels indicates the value obtained for the "whole slit" spectrum, (a) Profile of n e ; (b) top: profiles of T e ([N II]) (solid line) 
and T e ([0 in]) (dash-dotted line), bottom profile of the T e ([N ll])/T e ([0 III]) ratio; (c) profile of the dereddened flux ratio of [Fe III] 4658 
A; (d) top: observed flux of C II 4267 A (solid line) and O II 4649 A (dash-dotted line), botton: profile of the F(C II 4267 A)/F(0 II 
4649 A) ratio; (c) top: observed flux of [O III] 4959 A (solid line) and O II 4649 A (dash-dotted line), botton: profile of the F([0 III] 4959 
A)/F(0 II 4649 A) ratio; (f) top: 2+ abundances from RLs (solid line) and CELs (dash-dotted line), botton: profile of ADF(0 2+ ). All 
observed fluxes have units of erg cm -2 s~ 1 . Note that from panel a) to c) the profiles correspond to extractions of 1'.'2. From panels d) 
to f) the spatial increments are 4'.'8 wide. 



while the HG region and the area at the south of the po- 
sitional measurement ~95"of the slit position 2 are ionized 
by Her 36, the rest of the position 2 is ionized by another 
star, perhaps 9 Sgr. This hyp othesis is consistent wit h the 
results from previous work s of I Woodward et al l j 19861 ) and 
ISanchez fe PeimbertJ l|l99ll ). 

Finally, Figs[4jf)-[5ff) show the spatial variation of the 
2+ abundances derived from RLs and CELs, as well as the 
AD factor of this ion -ADF(0 2+ )- defined in its logarithmic 
form as 

ADF(0 2+ ) = Zo ff (0 2+ /H + ) RLs - /o ff (0 2+ /H + )cE Ls . (3) 

As in Figs[4jd)-(e) and[5jd)-(e), we present the spatial pro- 
file for the extractions 4'.'8 wide of the 2+ /H + ratio ob- 
tained from RLs and CELs. The ADF(0 2+ ) remains fairly 
constant along the slit positions 1 and 2 with average val- 
ues of 0.40±0.11 and 0.37±0.09 dex, respectively, which 
are in agreement with previous determinations obtained by 



lEsteban et all l|l999al ) -0.34 dex- and iGarcfa-Roias et al.l 
( 20071 ) -0.37 dex. The only relevant variation detected in 
the 2+ abundance profiles is at the HG region, where the 
2+ abundance decreases up to 7.6 dex due to the decrease 
of the ionization degree. 



5.2 M17 

Fig.[6]shows the spatial distributions of several nebular prop- 
erties along the slit position that covers two bright areas of 
M17 (see Fig. [1]). The density profile (Fig. (6^,) shows a peak 
between the positional measurements 30" and 90", a rather 
constant value between 90" and 180" and a decrease at the 
eastern edge. We have found a good agreement between our 
density det erminations and those at the positions 3 and 13 
studied by IPeimbert et~ai] l|l992h . which coincide with our 
positional measurements 60" and 125" (see iJSJ. In the case 
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Figure 5. Same as Fig.[4]fbr the slit position 2 of M8. The grey vertical band covers the apertures affected by stellar features from the 
emission of Her 36. 



of the temperature distributions (see Fig. [BJj) , we only have 
found variations of the order of our error bar. In the case 
of the T e (\Q ill]), we have fou nd values similar to those de- 
rived bv lPeimbert et al l (jl992r ) in their slit positions. On the 
other hand, the T,;([N n]) obtained by those authors shows a 
higher va lue at their pos i tion 3, which amounts to 11600 K. 
However, lEsteban et alj Jl999aT) also explored t he physical 
conditions at position 3 of iPeimbert et ah] (| 19921 ) and found 
a T e ([N n]) of about 8990 K, which is in agreement with 
our determinations within the errors. In Fig. [6jc), we can 
find the spatial profile of the [Fe ill] 4658 A line flux with 
respect to H/3, where we do not find any localized enhance- 
ment, which may be related to the presence of HH objects. 

The spatial distributions of the C II 4267 A and O II 
4649 A are presented in Fig. [BJd) , where we can notice an 
increase of the C II emission flux to the east of the positional 
measurement 140". This rise c an be related to a decr ease of 
the ionization degree. In fact, IPeimbert et al.l (|l992l ) find a 
change in the 2+ /0 + ratio between their slit positions 3 
and 13 of about 0.40 dex. The comparison of the [O ill] and 
O II emission line fluxes in Fig. H^e) does not show relevant 
features, and their ratio remains essentially constant along 
the slit. 

Finally, Fig. |6jf) shows the spatial profiles of the 2+ 
abundances derived from RLs and CELs, as well as their 



ratio, ADF(0 2+ ). As we can see, both abundance deter- 
minations show a similar decrease along the slit, while the 
ADF(0 2+ ) distribution present a constant value of about 
0.37±0.09 dex. This value is somewhat higher than that 
obtained in the "whole slit" spectrum, which amounts to 
0.31±0.07 dex, although the difference is consistent within 
the errors. 

5.3 NGC 7635 

The spatial distributions of nebular properties along the slit 
position of NGC 7635 are presented in Fig. [7] Each posi- 
tional measurement corresponds to an area of 3"x0'.'98. The 
results obtained for the integrated spectra covering sepa- 
rately the bright knots at the west of the central star, Kl 
and K2, and the rim of the bubble (see Fig.[T} are presented 
in 10 

In Fig. [7J a ) the density profile shows its highest value 
at the position o f knot Kl , wher e it reaches a value of 
about 2600 cm -3 . Rodriguez (1999) obtained the density in 
three zones of NGC 7635 and her positio n 1 coincides with 
our K l , as well a s the s lit position 3 of iTalent fe Dufourl 
(|l979T) . iRodn'guezl l|l999i ) derived a density value of about 
2800 cm" 3 , whic h is in agreement with our results, while 
ITalent fe Dufourl l| 19791 ) derived a density value of about 
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Figure 6. Same as Fig. [4] for the slit position of M17. The grey bands cover apertures where [Fe III] 4658 A (c) and O II 4649 A (d— f) 
were not detected. Positional measurements along the slit go from west to east (see Fig[Tj|. 



4200 cm . Moreover, iMoore et all (2002) also obtained a 
higher density values on the knots and the rim, but with 
larger uncertainties. Two local maxima can also be seen at 
the positional measurements corresponding to the knot K2 
and the rim. Between the knots and the rim, the density 
remains rather constant with an average value of about 100 
cm -3 and a standard deviation of about 20 cm -3 . 

In Fig. [TJb) , we can see that the electron temperatures 
do not show strong variations along the slit and are of the 
order of the error bar. The average values of both temper- 
atures amount to 8600 K for T e ([0 ill]) and 8100 K in the 
case of T e ([N n]) with standard deviations of about 390 K 
and 420 K, respectively. Their ratio shows an average value 
of 0.94 with a standard deviation of about 0.05. Our tem- 
perature d e termi n ations are in agreement with the results of 
iRodrfgued l|l999h . lTalent fc Dufourl (|l979T l and lMoore et all 
d2002h . 

Fig- El c ) snows the spatial profiles of the dereddened 
fluxes of [O ill] 5007 A and [O n] 7320 A with respect to H/3 
and normalized to their respective maximum emission. On 
the one hand, the [O n] spatial profile shows its maximun 
values at the position of the knot Kl decreasing progres- 
sively up to the positional measurement 110". On the other 
hand, the [O ill] distribution presents an inverse pattern 
with a progressive increase until the same positional mea- 



surement. Beyond this point, both spatial profiles do not 
show strong var iations. A similar r esult can be seen in the 
HST images of lMoore et al.1 (|2002h . The behaviour shown 
in Fig. [7J C ) is due- to a dramatic change of the ionization 
degree. The knots have a lower ionization degree than the 
gas beyond the positional measurement 110", which corre- 
sponds to the bubble. 



Finally, Fig. Etd) shows the heliocentric radial veloc- 
ity, Vhei , spatial distribution obtained from the Gaussian fit 
of the Ha line profile. We have obtained the most negative 
velocity at the position of the rim, which is about — 54±4 
an d in agreement with the mea n heliocentric velocity found 
by IChristopoulou et al.l l| 19951 ) for NGC 7635. This value 
is also similar to that of the molecula r cloud associated to 
the b ubble and the S162 complex (see IChristopoulou et al.l 
Il995l , and references therein). The expansion of the bubble 
is very slow, with velocities between 4 and 25 km s _1 and, 
therefore, our spectral resolution did not allow us to resolve 
the line splitting of the expanding bubble. In fact, the ve- 
locity profile along the slit shown in Fig. Etd) does not show 
any remarkable feature or velocity differences substantially 
larger than the uncertainties. 
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Figure 7. Spatial profiles of several nebular properties along the slit position of NGC 7635. Positional measurements along the slit go 
from west to east (see Fig[T}. The position of the knots, Kl and K2, and the RIM are indicated as well as the typical error of some 
variables. The horizontal dashed line in some panels gives the value obtained from the "whole slit" spectrum, (a) Profile of n e ; (b) 
top: profiles of T e ([N II]) (solid line) and T e ([0 III]) (dash-dotted line), bottom profile of T e ([N ll])/T e ([0 in]) ratio; (c) profile of the 
dereddened fluxes with respect to H/9 of [O III] 5007 A (solid line) and [O II] 7320 A (dash-dotted line) normalized to their respective 
maximum; (d) heliocentric velocity profile obtained from the centroid of Ho. 



Table 2. Physical conditions and 


ionic abundances 


1 of selected 


zones of NGC 7635. 










"Whole Slit" 




Kl 


K2 


RIM 


t 2 =0 i 2 >0 


n e ([S II]) (cm" 3 ) 


1900±150 


890±70 


180±30 


1080±90 


T e ([N II]) (K) 


8450±320 


8150±300 


8680±430 


8420±330 


T e ([0 in]) (K) 


8830±800 


8140±290 


8640±320 


8580±320 


c 2+b 


8.21±0.14 


8.46±0.08 


8.40±0.12 


8.44±0.12 


N+ 


7.55±0.04 


7.27±0.04 


6.87±0.05 


7.32±0.04 7.66±0.06 


o+ 


8.39±0.10 


8.28±0.10 


7.97±0.12 


8.24±0.10 8.45±0.17 


o 2 + 


7.32±0.10 


8.07±0.05 


8.08±0.05 


7.89±0.05 8.48±0.08 


Q 2+b 








8.48±0.09 


s+ 


6.37±0.04 


6.06±0.04 


5.53±0.04 


6.10±0.04 6.43±0.06 


s 2 + 


6.69±0.20 


6.90±0.09 


6.71±0.09 


6.77±0.09 7.40±0.12 


N+/O+ 


-0.84±0.11 


-1.02±0.11 


-1.10±0.13 


-0.91±0.11 -0.79±0.18 


o 2 +/o+ 


-1.07±0.15 


-0.21±0.11 


0.11±0.13 


-0.34±0.11 0.03±0.19 


S+/Q+ 


-2.02±0.11 


-2.23±0.11 


-2.44±0.13 


-2.14±0.11 -2.02±0.18 



a In units of 12+Zog(X +i /H+). 
b Determined from RLs. 
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6 CHEMICAL ABUNDANCES OF SELECTED 
ZONES OF NGC 7635 

In this section, we present the physical conditions and ionic 
and elemental abundances -derived from CELs and RLs- 
obtained from the integrated spectra that cover the knots 
-Kl and K2, the rim of the bubble as well as the "whole 
slit" spectrum. The results are shown in Tables [2] and [4] 

The physical conditions derived from the ratios of CELs 
are shown in Table [2] These values are in agreement with 
those obtained from the 3" extractions presented in Fig. 
as well as with the different determinations of the literature 
that we have cited in £15.31 We can notice that the tempera- 
tures calculated from the different indicators are very simi- 
lar within the errors, in agreement with the result obtained 
from the temperature spatial distribution along the slit po- 
sition (see Fig. [7p). The electron density and temperatures 
found for the "whole slit" spectra are basically averages of 
the knots and the rim values. 

We have derived ionic abundances of several ions (see 
Table [2J) • The + /H + ratio was obtained removing the tel- 
luric contamination from the [O n] 7320, 7330 A lines as it 
was explained in £|4. 21 A proper determination of the + 
abundances has allowed us to explore the strong variation 
of the ionization degree between the knots and the rim. As 
we can see in Tabled the ions with low ionization potential, 
such as N + , + and S + , show a gradual decrease in their 
abundance following the sequence: Kl, K2 and the rim. As 
it is expected, this decrease is accompanied by an increase 
of the abundances of ions of high ionization potential, how- 
ever the differences between K2 and the rim are not large 
in this case. In Fig. [8] we present the spatial distribution of 
several emission lines of ions with different ionization po- 
tentials: S+ (10.4 eV), H+ (13.6 eV) and 2+ (35.1 eV). 
The profiles are centered around the knots (Fig. [8^) and the 
rim (Fig. [SJa). As it is shown in Fig. [HI a), the [O ill] pro- 
file shows a completely different behaviour with respect to 
the other ions showing a progressive decrease from the posi- 
tional measurement 110" towards the west (left in Fig. [S£i). 
This behaviour is produced because most ionizing photons 
which are able to convert + to 2+ , are exhausted between 
the positional measurement 105" and 110" and, therefore, 
they do not get to penetrate the inner surface of Kl, which 
is facing the ionizing star. Then, Kl is a wall of dense ma- 
terial that maintains the nebula ionization-bounded at that 
precise location. In the case of K2, the spatial profiles of the 
three emission lines are rather similar, indicating that this 
feature is matter-bounded or conversely, it is a ionization 
front which surface is perpendicular to the line of sight. Fi- 
nally, the behaviour of the line profiles at the rim -whose 
surface is tangential to the line of sight- indicates that it is 
a matter-bounded feature. This is also consistent with the 
fact that T e ([0 ill]) and T e ([N n]) are so similar at the rim. 

It is important to note that we have detected and mea- 
sured the faint C II 4267 A RL (see Fig. [9]) in all the in- 
tegrated spectra, as well as five of the eight RLs of the 
multiplet 1 of O II (see Fig. [3JI) in the "whole slit" spec- 
trum, allowing us to derive C 2+ and 2+ abundances from 
RLs. Several authors ha ve previous l y dete rmined the C 2+ 
abundance in NGC 7635. iRodrfguezI |l999I) obtained values 
between 8.11 and 8.23 dex in her two slit positions which 
cover the group of knots at the south of Kl. On the bub- 
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Figure 8. Spatial profiles of several emission lines obtained from 
extractions with a size of 0'.'8 -the average seeing during the 
observations— around the knots (a) and the rim (b). 



ble rim, iMoore et all (2002) obtained a C 2+ abundance of 
8.63 dex. Our determinations lie in between those values. In 
addition, IMoore et all (2002) obtained the first determina- 
tion of the 2+ /H + ratio from RLs on the rim bubble, they 
found a value of 8.47 dex and a moderate ADF(0 2+ ) = 0.2 
dex. However, from our "whole slit" spectrum we obtain an 
ADF(0 2+ ) of about 0.59T0.10 dex, the highest value found 
so far in an H II region (see H7.ll) . 

In order to derive the total gaseous abundances of the 
different elements, we have corrected for the unseen ion- 
ization stages by using a set of ionization correction fac- 
tors (ICFs). Both adopted ICF values and total abundances 
are presented in Table [3] and Table [4l respectively. For 
carbon, we h a ve ad opted the photoionization models of 
iGarnett et al.l (|l998h to estimate the ICF(C + ). In the case 
of oxygen and nitrogen, we considered the classical assump- 
tions O = 0+ + 2+ and N+/N = Q+/Q respectively. Fi - 
nally, we have used the relation proposed bv lStasiriskal (|l97&f ) 
to derive the ICF(S 3+ ) and the S abundances. In general, 
the total abundances derived for Kl and K2 are rather simi- 
lar considering the uncertainties -except perhaps in the case 
of carbon- but higher t han those d e rived f or the spectrum of 
the ri m. Fo r the knots, Rodrigue3 (|l999h , iTalent fc Dufourl 
(|l9T9h and IMoore et al.l (|2002h determine a total O abun- 
dance of about 8.42, 8.72 and 8.78, respectively. Our O abun- 
dance determinations for Kl and K2 are only consistent with 
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Table 3. Set of ionization correction factors adopted in NGC 7635. 













"Whole Slit" 


Elements 


Unseen ion 


Kl 


K2 


RIM 


t 2 = 


t 2 > 


C 


C+ 


6.20±1.94 


2.03±0.55 


1.56±0.50 


2.34±0.66 


1.64±0.74 


N 


N 2 + 


1.09±0.33 


1.61±0.43 


2.28±0.73 


1.45±0.40 


2.08±0.93 


S 


S 3+ 


l.OOiO.Ol 


1.02±0.01 


1.07±0.04 


1.01±0.01 


1.05±0.05 



Table 4. Elemental abundances" in NGC 7635. 





C 


N 


O 


S 


Kl 
K2 
Rim 

"Whole Slit" l \ = ° 
t 2 > 


9.00±0.20 
8.77±0.14 
8.59±0.18 
8.80±0.17 
8.65±0.23 


7.58±0.14 
7.47±0.12 
7.23±0.15 
7.49±0.13 
7.98±0.20 


8.42±0.09 
8.49±0.06 
8.33±0.06 
8.40±0.08 
8.76±0.09 


6.86±0.14 
6.96±0.08 
6.77±0.09 
6.87±0.08 
7.46±0.11 



Abundances expected from abundance gradients 



Shaver et al. ( 1983") 


7.67±0.15 


8.60±0.15 


Afflerbach et al. C 1997) 


7.76±0.08 


8.56±0.11 6.67±0.07 


Deharvene et al. (2000) 




8.46±0.07 


Esteban et al. (2005)'' 


8.45±0.20 


8.64±0.12 


Cariei et al. (2005) c 


7.74±0.20 





Solar abundances 



Asplund et al. (2009) 8.43±0.05 7.83±0.05 8.69±0.05 7.12±0.03 

a In units of 12+log(X/H). 

b Determined from RLs. 

c Determined from CELs assuming t 2 > 0. 




\ (A) 

Figure 9. Sections of the integrated spectra around the emission 
line C II 4267 A for the knot Kl (top panel) and the rim (botton 
panel) of NGC 7635. 



those of iRodrfgue j l|l999l ). F or the rim, our O/H ra tio is 
somewhat lower than those of iTalent fe Dufourl (|l979l ) and 
iMoore et all (|2002h . both obtain a 12+Zop(0/H) = 8.49. The 
largest differences of abundances between the knots and the 
rim are found for carbon. The relative contribution of C + 
-and therefore the uncertainty of the total C/H ratio- is 
far larger in the case of Kl. One of the possible reasons of 
such difference may be the very different ionization degree 
of the knots and the rim, specially between Kl and the rim, 



where the 2+ /0 + ratio differs in more than 1 dex (see Ta- 
ble [5J. This strong change in the ionization degree, due to 
the ionization-bounded structure of the knots, involves an 
important variation in the ICF(C + ) (see Table[3]|. That fact 
can be affecting the suitability of the this ICF. 



7 DISCUSSION 

7.1 The abundance pattern of NGC 7635 

Assuming the validity of the tem p eratur e fluctua- 
tion hypothesis proposed by iPeimberU |l967) and that 
this phenomenon is related to the AD problem (see 
iGarcfa-Roias fc Este ban 2007!), we have estimated a t 2 pa- 
rameter from the ADF(0 2+ ) found for the "whole slit" spec- 



tra, fo llowing the formalism outl ined by Peimbert fc Costerol 
l| 19691 ) and equations (8)-(ll) of lPeimbert et al.l l|2004l ). The 
value of the t 2 parameter representative of the 2+ zone is 
0.071±0.009. We have adopted this t 2 parameter in order to 
calculate the ionic abundances, the same set of ICFs (see 3St 
and the total abundances for the "whole slit" spectrum un- 
der the presence of temperature fluctuations. These values 
are also included in Tables O and U 

The "whole slit" abundances can be considered as the 
representative average values for NGC 7635 and, therefore, 
it seems appropriate to compare them with those expected 
from the Galactic abundance gradients. We have included 
in Table [4] the expected abundances obtained from some of 
the gradients determinations available in the literature, as 
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well as the solar values. U sing the dist a nce t o NGC 7635 
of 2.4±0.2 kpc adopted bv lMoore et~ai] \200j ) and assum- 
ing a galac tocentric distance of the Sun of 8.0±0.5 kpc 
jReidll 19931 ). we estimate a galactocentric distance of about 
9.2±0.5 kpc for NGC 7635. As we can see in Table H the 
comparison of the total abundances for both possible val- 
ues of t 2 and the abundance gradients is rather puzzling. 
Firstly, the C/H ratio obtained from the C n 4267 A RL 
is larger than expected by the Galactic abundance gradi- 
ent of this element de rived from observations of H II regions 
jEsteban et all 120051 ) for any of the two values of the t 
parameter considered in the table. The difference between 
our C abundance determinations and the expectations of 
the abundance gradients is specially large in the case or the 
knots, but marginally consistent within the errors in the 
case of the rim. This could suggest that -as we outlined in 
Section [6]- this disagreement may be due to the unsuitabil- 
ity of the ICF(C + ) used, at least in the case of the knots. 
In the case of the "whole slit" spectrum, the C 2+ /0 2+ ra- 
tio with respect to the excitation degree -0 2+ /0 + ratio- of 
the ionized gas is similar to the values reported for other 
Galactic H n regions, as S 311 (|Garci'a-Roias et al.l 12005? ) 
and, therefore, the ICF(C + ) seems to be not the reason of 
the large C/H ratio in this particular case. The expected 
values of the Galactic O/H and N/H gradients are just in 
between our determinations for t 2 =0 and t 2 > 0. In partic- 
ul ar, the O abun d ance e xpected in the abundance gradient 
of lEsteban eTakl (|2005h -based on the flux of O II RLs- 
is marginally consistent with our determinations for t 2 > 
within the errors. Finally, the S/H ratio calculated for any 
value of the t 2 paramet er is always large r than that expected 
from the S gradient of lAfnerbach et all l|l997t) . 

The lack of consistency between our abundance deter- 
minations and the Galactic abundance gradients expecta- 
tions for NGC 7635 suggests that perhaps the standard 
methods for deriving physical conditions and chemical abun- 
dances are not giving the correct values for this object. The 
problem seems to affect the spectra of all selected areas: 
the knots, the rim and the "whole slit" ones. Firstly, the 
carbon abundance is too high. It is very unlikely that the 
nebula has suffered a pollution of C processed by the ioniz- 
ing central star because it is still in the main sequence. On 
the other hand, due to the fact that Kl and K2 also show 
a very intense C n 4267 A RL and high C 2+ /H + ratios, 
the hypothetical chemical pollution event should also affect 
these structures that belo ng to the larger emission complex 
S162 (|Moore et all 120021 ) . which is outside the bubble of 
NGC 7635. Another possibility of a local increase of C may 
be due to the destruction of carbon-rich dust by the shock 
associated with the expanding windbl own bubble. A pos - 
sibility that was alre ady suggested by I Moore et all (2002). 
lEsteban et al.l (| 19981 ) estimated that the destruction of all 
the carbon locked-up onto dust grains in the Orion Nebula 
would increase the measured gas-phase C/H ratio in only 
about 0.1 dex. An increase of such magnitude would allevi- 
ate somehow the difference between the observed values and 
those expected by the abundance gradients, but does not ac- 
count for it completely. All these problems with the flux of 
the C II 4267 A RL and the high C/H ratio, specially at 
Kl, makes us to consider that the emission flux of that line 
is enhanced by an unknown mechanism. We speculate that 
perhaps this mechanism could be related with the strong 



ph otoevaporative flow from the surfaces of the knots found 
bv lMoore et"afl (|2002h from HST images of NGC 7635. For 
these authors, structures as Kl are the ionized edge of a 
mass of neutral material associated to the S162 complex, 
which is ionized by the central Of-type star. The photoe- 
vaporative flow seem to be interacting with the stellar wind 
as it is suggested by the presence of a n emission loop be - 
tween the knots and the central star (|Moore et al ] |2002h . 
Further deeper and higher spectral resolution spectroscopi- 
cal observations of these features would be needed to shed 
light onto this problem. 



7.2 Comparison with the Orion Nebula 

The Orion Nebula is the nearest H II region and there- 
fore the most suitable object to explore the spatial vari- 
ations of the nebular properties as well as the abun- 
dance discrepancy problem at the highest spatial resolution. 
M esa-Delgado et all l|2008h obtained the spatial distribu- 
tions of several quantities over five slit positions in the Orion 
Nebula with an angular resolution set to l'/2. They found 
significant small-spatial scale variations of density, temper- 
ature as well as the ADF(0 2+ ), most of them related to the 
presence of morphological structures such as HH objects, 
proplyds and ionization fronts as the Orion bar. 

In this study, we have tried to obtain a similar dataset 
for other bright -but more distant- Galactic H II regions. 
However, the inventories of proplyds and HH objects in M8 
and M17 are far less complete than for the Orion Nebula. Be- 
cause of its different nature, there is no detection of these ob- 
jects in NGC 7635. Discovering proplyds in other H II regions 
-apart from the Orion Nebula- is a difficult task considering 
their larger distances and the problems related to the dis- 
crimation between true proplyds and fragmen ted portions 
of molecular clouds fsee lDe Marco et al.l 120060 . The single 
ca ndidates to prop l yds in M8 and M17 have been reported 
by IStecklum et alj Jl99Sf) -marked as UC in Fig. [f]- and 
|Pe Marco et all l|2006h ~ analyzing narrow-band images taken 
with HST, respectively. On the other hand, HH objects are 
direct manifestations of the interaction of gas ejected by 
a young star with its surroundings. The existence of HH 
objects in M8 was firstly mentioned by iReipurthl l|l98lf ). 
who discovered HH 213 from images with the 1-m telescope 
class at Las Campanas Observatory. More recently, the cat- 
alog of HH objects in M8 has i ncreased with the spectro- 
scopic confirmation of HH 87 (lArias et al J l2006h and the 
new five outflows reported by iBarba fc Ariasl |2007l ) using 
optical narrow-band imaging with the Wide Field Imager at 
the 2.2-m telescope at La Silla Observatory. In the case of 
M17, there is not available literature about the presence of 
HH obje cts in this H H region, e xcluding a possible microjet 
found bv lDe Marco et all (|2006l 1. 

Assuming a distanc e to the Orion Nebula of 436 pc 
l |0 'Dell fc Hennevl [2008 ) , the linear sizes of the proplyds 
and HH objects located at the Huygens region lie between 



0.003 and 0.007 pc and 0.006 and 0.02 pc, re s pectiv ely. The 
extractions analyzed by iMesa-Delgado et al.l (120081 ) had an 
angular size of l'.'2, that corresp onds to a linear resolution o f 
0.0025 pc. Given that resolution. IMesa-Delgado et all (|2008T ) 
could obtain enough spatial sampling to detect the spatial 
variations of the nebular properties associated to those mor- 
phological structures. However, the H II regions we analyse 
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in this paper are at larger dista nces than the Orion Nebula: 
1.25 kpc (M8 1 Arias et alj|2006h. 16 kpc (M17|Povich et all 
\200H ) and 2.4 kpc (NGC 7635 iMoore et alJl2002h . At those 
distances and with the size of the extractions used for each 
H II region in this work, the minimum linear size we can 
resolve is 0.0072 pc, 0.0093 pc and 0.035 pc in the case of 
M8, M17 and NGC 7635, respectively. As we can see, the 
linear resolutions that we achieve in our spatial analysis are 
a factor 3, 4 and 14, respectively, lower than in the case of 
the Orion Nebula. It is clear that spatial variations of the 
physical and chemical properties associated with proplyds 
with linear sizes similar to those of the Orion Nebula are 
very unlikely to be detected with our observations, except 
perhaps in the case of M8, where very large proplyds are of 
the order of our small extractions. With these comparatively 
large resolution elements, the emission of the proplyds -and 
any associated variation of the local properties of the ion- 
ized gas, if they exist- would be diluted with the emission 
of the ambient nebular gas. In the case of HH objects, struc- 
tures similar to those observed in the Orion Nebula would 
be resolved with the size of the extractions used in the cases 
of M8 and M17. For example, knot B of HH 870 in M8 has 
a diameter of about 0.015 pc, similar to those found in the 
Orion Nebula. 

If those object are present in all H II regions, we have 
estimated that the spatial resolution needed to resolve pro- 
plyds similar to those of the Orion Nebula should be about 
0'.'42 for M8, 0733 for M17 and 0'.'22 in the case of NGC 7635, 
which are of the order or smaller than the best seeing at- 
tainable with ground-based telescopes. 



7.3 A possible candidate to Herbig-Haro object in 
M8 



In §571] we have presented several evidences that suggest 
the presence of a new HH object in M8. It is located ap- 
proximately 16" east and 44" north from Her 36 and can 
be seen as a diffuse relatively bright arc near an ioniza- 
tion front (HH? in Fig. [T]). This candidate to HH object 
presents a density 1.4 times larger than the adjacent back- 
ground, somewhat lower than the density contrast found for 
the much prominent HH 870, which amounts to 2.3. 

As we mentioned in H5.ll the HH objects typically show 
a strong emission of lines emitted by low ionization poten- 
tial ions such as Fe 2+ , S + or + . The spatial distribution 
of [Fe ill] 4658 A shown in Fig. gjc) indicates that this line 
is enhanced by a factor of 1.3 with respect to the surround- 
ing background gas. This increase of the dereddened flux of 
[Fe ill] 4658 A is also lower than that shown by HH 870. The 
[S n]/Ha ratio is a good indicator t o discrimina t e betw een 
shock-excited and photoionized gas. lArias et al. | (|2006h ob- 
tained a map of [S ii]/Hq ratio of the central part of M8 
nebula, where they reported high values of this indicator at 
the three nebular knots associated with HH 870. Using our 
data, we find that the [S ll]/Ha ratio increases by a factor of 
1.8 at the position of HH 870. In the case of our candidate 
to HH object, we have obtained an increment o f about 1.3, 
simila r to that we can estimate from the map ofH rias et al.l 
(2006). In fact the shape of the possible HH o bject can be 
well id entified in the [S n] and [S ii]/Hq maps of lArias et al.l 
(2006), while it cannot be noticed in their Ha map. This 
behaviour is different to that we can observe in the ioniza- 



tion front near the candidate, which can be clearly seen in 
all maps. Further spectroscopic observations at higher spec- 
tral resolution would be necessary to study the kinematics 
of this object in order to confirm its true nature. 



8 CONCLUSIONS 

In this article, we have carried out long-slit spectrophotom- 
etry at intermediate spectral resolution of the Galactic H II 
regions M8, M17 and NGC 7635. The one-dimensional spec- 
tra were extracted with a resolution of l'.'2 for M8 and M17, 
and 3" in the case of NGC 7635. Additional extractions with 
a spatial size of 4" 8 were necessary in order to measure the 
faint C II and O II RLs in M8. We have studied the spatial 
distributions of a large number of nebular quantities along 
several slit positions covering different morphological struc- 
tures such as HH objects, ionization fronts or bright knots. 
The studied quantities were c(Hj9), n e , T e ([0 in]), T e ( [N il]), 
the observed and dereddened flux of several emission lines 
([Fe m] 4658 A, C n 4267 A, O n 4649 A, [O in] 4959,5007 
A and [O n] 7320 A) and the 2+ abundances derived from 
CELs and RLs, as well as the difference of both determina- 
tions, ADF(0 2+ ). 

The density spatial distributions show a large range of 
variation across the different slit positions. We have found 
local maxima associated with the HG region, HH 870, knots 
and regions with a high surface brightness. On the other 
hand, the temperature spatial profiles do not show impor- 
tant variations related to the cited structures. The temper- 
atures obtained from the different indicators present the 
classical behaviour for M8 and M17: those derived from 
[N n] lines are higher than those derived from [O ill] lines, 
as expected for ionization-bounded nebulae. In the case of 
NGC 7635, both temperatures seem to be very similar con- 
sidering the error due to the own structure of the Bubble 
nebula which is matter-bounded. We have also explored the 
spatial behaviour of the ADF(0 2+ ) along the slit positions 
of M8 and M17, which remains rather constant, finding val- 
ues between 0.3 and 0.5 dex with an average error of about 
0.1 dex. 

We have analysed the physical conditions and chemi- 
cal composition of four additional extractions of NGC 7635: 
three of them centered on the knots -Kl and K2- and the 
rim of the bubble; and the last one corresponding to the 
"whole slit" spectrum. On the one hand, the comparison of 
the 2+ /H + ratio determined from CELs and RLs in the 
"whole slit" spectrum produces an ADF(0 2+ ) of about 0.59 
dex. Assuming that the AD problem is related to tempera- 
ture fluctuations, we have obtained a t 2 parameter of 0.071. 
We have found a puzzling pattern in the total abundances 
derived for NGC 7635. The total abundances obtained for 
the knots are slightly higher than those of the rim. The to- 
tal abundances were compared with those expected by the 
Galactic abundance gradients, finding that there are discrep- 
ancies, specially in the case of C. We suspect that C II 4267 
A RL may be abnormally enhanced in NGC 7635 due to 
an unknown physical process, whose investigation is outside 
the scope of this paper but deserves further more detailed 
observations. 

Comparing o ur ob servations with those of 
M esa- Delgado et al] ((2008), we conclude that proplyds 
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-and the associated variations of the local properties of the 
gas- with linear sizes similar to those found in the Orion 
Nebula cannot be resolved with the observations reported 
in this paper. Angular resolutions of the order or smaller 
than the minimun seeing reached from the ground-based 
telescope would be needed to distinguish the presence of 
proplyds. 

Finally, we have found several evidences that point out 
to a possible new candidate to HH object in M8. This 
new object is located 16" east and 44" north from Her 36 
where we have found enhancements in the spati al profile of 
[Fe in ] 4658 A and in [S ii]/Hq map presented bv lArias et al.l 
l|2006l i. 
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(Hi) We have not included the errors associated with H/3 
observed fluxes which remain between 3% and 5% among 
the different slit positions. 

(iv) We have only listed [O ill] 4959 A and [N n] 6548 A neb- 
ular lines. [O in] 5007 A and [N n] 6583 A can be obtained 
from the first ones using their theoretical relation: 

J([0 in] 5007)//([O m] 4959) = 2.88 (Al) 

and 

/([N n] 6583)/J([N n] 6548) = 2.92. (A2) 

(v) In the cases of M8 POS1 and POS2, it should be re- 
minded that additional extractions of 4'.' 8 wide were ob- 
tained to achieve proper flux measurements of O II RLs (see 
£[3j . In this sense, an aperture of 4'.'8 wide represents the av- 
erage value of four apertures extracted with an angular size 
of l'.'2. The vertical lines in the left-hand side of the O II 
RLs cover the four apertures of l'.'2 wide. 



APPENDIX A: EMISSION LINE FLUXES 

In Table Al, we present the aperture number (Column 

1) , H/3 observed fluxes in units of erg cm 2 s _1 (Column 

2) , dereddened flux line ratios of the main emission lines 
(Columns 3-14) per slit position in untis of /(H/3) = 100 
and the average extinction coeficient (last row of each slit 
position). Several notes should be considered in order to un- 
derstand Table Al: 

(i) The aperture number (Ap) is the identification num- 
ber of each one-dimensional spectra extracted and it is 
related to the positional measurement in arcseconds as 
Api n i t i a i-|-Ap 3 i Ze xAp. Apiniuai corresponds to the position 
in arcseconds after the edges of the CCD were discarded 
for the extraction of one-dimensional spectra (see extrac- 
tion procedure in Sj3]). Ap ini ti a i amounts to 34" in M8 POS1 
and POS2, 10" in M17 and 82" in NGC 7635. Ap size is the 
extraction size in arcseconds along the spatial direction (see 
the selected sizes in Sj3]). 

(ii) In order to interpret correctly the Ap column in M8 
POS2 we should consider that apertures 65, 66, 67 and 68 
were removed in this position due to contamination by stel- 
lar emission (see fj3j> . 



